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Abstract 

The discrete Painleve III equation is investigated based on the bihnear formahsm. 
It is shown that it admits the solutions expressed by the Casorati determinant whose 
entries are given by the discrete Bessel function. Moreover, based on the observation 
that these discrete Bessel functions are transformed to the g-Bessel functions by a 
simple variable transformation, we present a g-difference analogue of the Painleve III 
equation. 



1 Introduction 



Recently, several interesting ordinary difference equations have been proposed as discrete 
analogues of the Painleve equations. Among them, the discrete Painleve I(dPi) and Il(dPn) 
are naturally derived in calculating the partition functions of some models in the theory of 
two-dimensional quantum gravity In other words, they have close relationship with 
the theory of orthogonal polynomials. This example implies the significance of discrete 
Painleve equations. Moreover, they are expected to possess rich mathematical structures. 
Besides this derivation, the discrete Painleve equations can be obtained in several con- 
textii|,0,i. 

It is known that the discrete Painleve equations admit several properties analogous to 
continuous ones, such as the existence of Lax pair, Backlund transformation, coalescence 
cascade and singularity confinement (SC) property [p|,|6|,[7|,p|, p^ , pT| , p!2| . However, are these 
properties sufficient to regard the equations really as the discrete Painleve equations? 

In the continuous case, the Painleve equations are "integrated" in the sense that they 
are reduced to linear integral equations by means of inverse scattering transform from their 
Lax pairs|l^. In the discrete case, however, the dPi equation is the only known example 
which is "integrated" based on its Lax pair. 

In the continuous case, series of exact solutions which are expressed by determinants 



are obtained by using the Backlund transformation |jTj]. However, no attempt has been 
done to construct such solutions for the discrete case. 

The coalescence cascade is a property of the Painleve equations which is related to 
the confluence of singularities of linear ordinary differential equations in the theory of 
monodromy preserving deformation. Hence, this property is considered to have no direct 
relationship to the structure of their solution spaces. 

The SC is a concept which was extracted as a common property from the class of 
integrable ordinary difference equations proposed by Quispel et al.|TH] In this context, "in- 
tegrable" means the possession of conserved quantities. Grammaticos et al. have employed 
this property as an integrability criterion of difference equations^]. By applying this test 
to non-autonomous equations, they have proposed several discrete Painleve equations^. 
Since the SC might be regarded as the discrete analogue of the Painleve property, it pro- 
vides us with a useful and powerful tool to detect "good" equations. However, we do not 
know whether the equations which pass the SC are really "integrable" in the above sense. 

The Painleve equations are characterized by the property that they do not possess 
movable singularities except for poles. It is also a significant property that they admit series 
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of exact solutions expressed by determinants for special values of parameters, although their 
solutions are transcendental in general. We believe that the existence of such solutions is 
a direct reflection of the "integrability" of equations. From this viewpoint, we consider 
that the preservation of the structure of solution space is the most important for the 
discretization of Painleve equations. 

In the previous paper, we have shown that the dPn admits a class of solutions expressed 
by the Casorati determinants of discrete Airy function They reduce to the Wronskian 
solutions of Pii in the continuous limit. In this paper, we investigate the solutions of special 
function type for the discrete Painleve III (dPni) equation. We will show that they are 
expressed in terms of the Casorati determinants of the discrete Bessel functions. The dPm 
occupies a particular place among the discrete Painleve equations, since this is the simplest 
equation which is derived only through the SC. Hence, the existence of such solutions may 
be a good check of the validity of SC. 

Moreover, we will show that the discrete Bessel functions in our result are transformed 
to the g-Bessel functions by a simple independent variable transformation. Based on this 
observation, we present a g-difference analogue of the Pni. 

In section 2, we present the explicit form of exact solutions for the dPni- The derivation 
is discussed in detail by means of the bilinear formalism in section 3. In section 4, we present 
a g-difference analogue of the Pm. Finally in section 5, we give concluding remarks. 



2 Solutions of dPni 



In this section, we investigate the solutions of dPin[§], 

Wn+lWn-l = ■ , (1) 

where a, /3, 7, (5, and A are parameters. Before going to the discrete case, let us remind of 
the special function type solutions for the Pm equation. 



d?u 1 / du\ 1 1 , n ^, 5 

= - — ---r + +I3) + W + - , 

dx"^ u \dx X ax X u 



(2) 



where a, /3, 7 and 5 are parameters. It is known that eq.(|D admits the solutions expressed 
by the r function ||14||, 



r; 



N 





D J, ■ 






D J, 


D'J, ■ 


■ 


J. 








'J. 



(3) 
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where D = x-^ and is the Bessel function of degree i^. The dependent variable u in 
eq.(ll) is expressed as 

„,.f.ogirU^. (4) 



Equation (^) gives the particular solutions of eq.(0) with the parameters, 

a = 2{N-u), P = 2{u + N + 1), 7 = 1, S = -1 . (5) 
To obtain the particular solutions of dPm eq.(|I]), we start from the discrete Riccati 



equation[17|. In this case, it is easily shown that if Wn satisfies 



aWn + A" 

Wn + d 

where a and d are constants, then Wn gives a solution of eq.(|l]). Assuming the fractional 
form, 

Wn= , (7) 

and substituting this into eq.(|), we get 

Wn = a , (8) 

Gn+2 + (a - d)Gn+i + (A" - ad)Gn = . (9) 

On the analogy of the continuous case, it is expected that eq.(D is a discrete version of 
the Bessel equation. However, it is not trivial to identify eq.(^ with the Bessel equation. 
Instead of eqs.(|^) and (^), we consider 

Jy[n) 

Mn + 2) - (g^ + q-^Mn + 1) + {1 + (1 - g) Vn^.H = ■ (11) 

We can show that if Ju{n) satisfies eq. ([Tl|) , then Wn in eq.(p!0D gives the solution of dPm 
eq.(|l]) with the parameters, 

a = -l, p = [q- - q—^){l - q)\ 7 = ^^^, 5 = q^ - q-^ , \ = q' . (12) 

Continuous limit of eq.(0) is given by putting g = l + e, ?t,= |- and taking the limit of 
e — ^> 0, which yields 

^ + (e--.^)J. = 0. (13) 
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Equation (|T3D reduces to the Bessel equation in ordinary form if we put x = e^. In this 
sense, eq. (pTf) is regarded as the discrete Bessel equation. We note that the same results 



for this simplest solution is also derived in ref. . 

The critical point is that the above result can be extended to the N x N determinant 
solution. We here give only the results, leaving the derivation in the next section. We 
consider the r function. 



Jv{n) 



Ju{n + 3) 



J^{n + 2N -2) J^{n + 2N -I) 
where Ju{n) satisfies the following contiguity relations. 



J,{n + N-l) 
J,{n + N + 1) 

J^{n + 3N-3) 



(14) 



Mn) - q-'Mn + 1) = -(1 - g)g" J,+i(n) 
Mn) - q^Mn + 1) = (1 - g)g"J,_i(n) . 



(15) 
(16) 



We notice that the discrete Bessel equation (|TT]) is derived /^from above contiguity relations. 
Then we can show that 



- q 



rN 



" rj^+i(n)rr^(n + l) 
gives the solutions of the dPm, 



(17) 



Wn+lWn-l 



w'i + 5wn + a 



with the parameters. 



a 



AN 

-q , 



P = {q 



v+N 



-u-N-2\ ^8N 



q^^il-qf 



(19) 



It is remarkable that the r function has a similar asymmetric structure to that of dPn|16 



Namely, the index of elements of the determinant, n, increases by two in one direction 
while by one in another direction. 



3 Derivation of Results 

In this section, we show that the r function (|1^) actually gives the solution of eq. (0) with 
the parameters ([I9|), through the dependent variable transformation ([T7|). First, we give 
the bilinear difference equations satisfied by the r function (p^if) . Secondly, we show that 



eq.(18) follows from the bilinear difference equations. 
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3.1 Bilinear Difference Equations 

The bilinear difference equations for tlie r function (|n|) are given as 

rN+iin) r'^,^\n + 1) - q'^-"" r'^^^,{n + 1) r^^+^H 

= -(l-g)g"+2^r;^+\(n)rj^(n + l), (20) 

rp-_,\{n) r'^,{n + 1) - g'^-^^ (n + 1) rj^(n) 

= (l-g)g"+2^r;^^,(n)r;^+i(n + l), (21) 

r^N^\n + 3) - g-'^-^r]^+i(n + 1) r^^^\n + 2) 

= -(l-g)g'%+\(n)r;;(n + 3), (22) 

r^+iH + 3) - g''-^+V;^+\(n + 1) r^^{n + 2) 

= (l-?)«+iWrj^^'(^ + 3), (23) 

r]^^,(n + 2) r;^+i(n + 1) - ^^^(g'^-^ + g-'^+^)r;^+i(n + 1) r^^^\n + 2) 

+ g^^{l + (1 - g)^g^"}r;^+iH ^^^^'(^ + 3) = . (24) 

In tlie following, we give a detail of the proof that eqs.(^)-(|2^) hold for the r function 

(0) - 

(1) Difference Formulas 

The Pliicker relations are the quadratic identities for determinants some of whose 
columns are shifted. Bilinear equations are obtained from the relations. Therefore the 
first step is to construct several formulas which relate the "shifted determinants". For 
usual integrable systems, such as the KP, Toda lattices and discrete KP, the entries of 
their determinant representation of r functions satisfy quite simple linear differential or 
difference equations, and hence we can easily construct the differential or difference formu- 
las for the r functions (see, for example, |lT8| , p!9| , pO| , ^ ) . On the other hand, in this case. 



the entries of the r functions satisfy linear, but variable-coefficient difference equations as 
shown below, which makes the calculations somewhat technical. 

We now consider the r function ( p!^ provided that the element Ju{n) satisfies the 
contiguity relations (|T3p and (0). Let us introduce a notation for simplicity, 

r;^(n) = |0„l„---,iV-2„iV-l,| , (25) 
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where k^, stands for the column which begins with J^{n + k), 

( J^{n + k) \ 
J^{n + k + 2) 
J^{n + k + 4) 



V 



(26) 



Here the height of the column k^ is A^. However in the following, we use the same symbol 
ki, for the bigger determinant T^^i{n) because there is no possibility of confusion. So the 
size of ki, should be read appropriately case by case. By adding (z + l)-th column multiplied 
by —q'"^ to i-th column and using eq. ([T5|) for i = 1 ~ (A^ — 1), we have 



«r -, (2n + ]V-2)(JV-l) 

Ju+iin) 



X 



gV^+i(n + 2) 



q^^-^Ju+i{n + 2N-2) ■■ 

..N-l (2"+3iV-2)(iV-l) 

{-(1-?)} Q 



J,+i{n + N-2) J,{n + N-l) 

q^J^+i{n + N) jJ{n + N + l) 

g2^-2j^+i(n + 3Ar-4) J^(n + SAT - 3) 



(27) 



X 



Ju+i{n) 
J^+i(n + 2) 

J^+i(r2 + 2A^-2) 



J,+i{n + N-2) J,{n + N-l) 

Jy+iin + N) q-^J^{n + N +1) 

J^+i (n + 3 - 4) g-27v+2 + 3 AT - 3) 



,r -, (2n + 3]V-2)(iV-l) 

{-(1 - g)}^-ig' |0,+i, ■ ■ ■ , A^ - 2,+i, A^ - i;. 



where A;^ stands for the column, 



k' 



( Ju{n + k) \ 
q-^J^{n + k + 2) 
q-^J,,{n + k + A) 



\ 



i2t 



J 



Hence we obtain a difference formula. 



«r , -, (2n+3JV-2){]V-l) , 

|0., • • • , AT - 2., AT - = {-(1 - g)}-^+^g-^ T^^-\n). (29) 

Let us construct another difference formula. In eq.(^), multiplying A^-th column by 
q'", adding (A^ — l)-th column multiplied by —(1 — g)g"+^~^ to A^-th column and using 
eq.(|l5D, we have 



q "r^in) = {-(1 - g)}'^ 'q 
Ju+i{n) 
q^Ju+i{n + 2) 



X 



g'^-V.+i(n + 2Ar-2) 



J^+i(n + A^-2) J^{n + N-2) 

q^Ju+i{n + N) Juin + N) 

q^^-^Jy+iin + m -A) J^{n + ^N -A) 
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(2n+3iV-2)(iV-l) 
2 



X 



Ju+i{n) 
J,+i{n + 2N-2) 

{-(1-g)} ^ ^ 



J,+i{n + N-2) 
J,+i{n + N) 



J,{n + N -2) 
q-^Juin + N) 



J^+i (n + 3 - 4) + 3iv - 4) 

.,iV-2,+i,iV-2'J. 



Hence we get a difference formula, 

|0,, • • • , iV - 2,, iV - 2:,_i| = {-[l-q)y-'^q 
Similar calculations give the following formulas by using eq. ([l6|) instead of (|l5l), 



-jV+1 - (^"+3iV-^)(^-l) -^+l 



|0,,---,iV-2,,iV-i; 



|0. 



iV-2,,iV-2',+J 



.J , , (2n+3]V-2)(JV-l) , , -, , -, 



(30) 

(31) 
(32) 



So far we have constructed the difference formulas of r for "horizontal shifts" . Next we 
consider formulas for "vertical shifts" (shifts of rows in eq.(|T^). For notational convenience, 
we take the transpose of r. 



r;^(n) 



Ju{n) 
Ju{n + 1) 



Mn + 2) 
Mn + 3) 



J,{n + N-1) J,{n + N + l) 
|0„2,,---,2iV^4„2iV^2,| . 



J^{n + 2N -2) 
jJin + 2N -I) 

J^(n + 3Ar-3) 



(33) 



Here ku stands for the column, 



/ Jy{n + k) \ 
J^{n + k + l) 
J^{n + k + 2) 



V 



(34) 



Using the contiguity relations (|I5|) and (0), we can verify 
2.-g^'^(l + (l-g)V") 6. 

g-H(l + q-"'^')J.+iin + 2) - (1 - g^)^— i J,^,(n + 1)} 



(1 - q)q 



g-2{(l + q-"'+')J.+iin + 3) - (1 - q^)q-^-' J,+,in + 2)} 



v 



■ (35) 



In eq.(^), we first add j-th column multiplied by —q^'^{l + (1 — g)2g2n+4j 4^ _|_ ]^^_^]-^ 
column for j = (A^ — 1) ~ 1 and use eq.(^5|). Secondly, we multiply i-th row by g*^^/ (1 + 
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^_2^+2{i-i)) for i = 1 ~ A^. Finally adding i-th row multiplied by (l-g2i)g-'--i/(i+g-2^+2i) 
to [i + l)-th row for i = 1 ~ (A^ — 1), we obtain 



{{l-q)q 



n+3u-l\-N+l 



y 



N-1 



1 



(36) 



where 



k,, 



Ky{n + k)i 



\ 



(37) 



and Ky{n)i is defined by the following recursion relation, 



KJn) 



KJn)o 



1 



1 



^{q^J,{n + + (1 - q^')q-'-'K,in),^,), for z > 1 



-2u 



JJn) 



{31 



Let us derive another difference formula. We multiply 1st column of T^{n) in eq.(|33|) by 
— g^'^(l + (l— g)^g^") and add 2nd column to 1st column. Then we follow the same procedure 
as above calculation, that is, add j-th column multiplied by —q^^{l + (1 — q)'^q^'^^'^^~^) to 
(j + l)-th column for j = (A^ — 1) ~ 2, multiply i-th row by q^^^/{l + g"2i/+2(j-i)^ i^^^ 
z = 1 ~ iV, and add z-th row multiplied by (1 — q^'')q~^^^ / {I + g-2v+2j^ _|_ iy\^\^ j~q^ 
for z = 1 ~ (A^ — 1). Then we obtain 



llj'+i) 2^, 3i^+i, ■ ■ ■ , 2A^ — 3i^+i| 
{(l-g)g"+3'^-i}-^+ig-^^ 



iV-l 



n 1 1 „-2i^+2fc 

k=Q J- -t- (/ 



X (V'^){l + (l-g)V"} rj^(n). 



(39) 



We have constructed six difference formulas (^), (^), (0), (g^), (^ and (^. In 
the following, using these difference formulas for r function, we prove that the Pliicker 
relations reduce to the bilinear difference equations (^)-(^). 
(2) Bilinear Identities 
A. Equations and (0) 

First we derive eq. (pOD . Consider the identity of determinant, 

0, I, ■■■ N-l, \ I 01 











(40) 



8 



where fcj, and denote columns defined in eqs.(2C) and (pS]), means an empty matrix 
of relevant size and 

■ U+1. (41) 
Applying the Laplace expansion on the right-hand side, we have 



= |0,,---,iV-l,,iV:_i| X |l„---,iV-l„iV„0i| 

- \1„---,N -U,N„Nl_^\ X |O„---,iV-l„0i| 

- |0„ ■ ■ ■ , - 1„ A^.l X |1„ ■ ■ ■ , - 1„ Nl__^, 0i| 



(42) 



We note that eg . (14^) is nothing but one of the Pliicker relations. This identity is rewritten 
in terms of r by using eqs. (pS]), ( pO]) and (|30|) as 



-(l-g)g"+^^r;^+i(n)r;^-^(n+l) 



(43) 



which is essentially the same as eq.(^). Equation (^) is derived from the same identity as 
eq.(^2|) except that Nl_^ is replaced by A^^+i and that eqs.(PTD and ( P^ are used instead 



of eqs. (H) and (0). 

B. Equations (|2|) and ([23D 

We next derive eq.(^). Replacing 0i in eq.(^2D by 







Vo/ 



we have 



(44) 





■,N- 


1., 


A^' 


-1 ^ ll') 


■■■,N 


- l.,Ar„< 


P2\ 




■,N- 












P2\ 


- |o.,-- 


■,N- 


1., 




X I'^u, ■ ■ 


■,N- 




P2\ 



(45) 

which can be rewritten in terms of r again. Here we give notes on the treatise of the deter- 
minant including the column 02- As an example, let us consider ■ ■ ■ , N — 1^, Nl_^, 02 1- 
Expanding by (A^ + l)-th column, we have 



|l„---,A^-l„Ar^_i, 
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Juin + l) 
J.(n + 3) 



J,{n + N-l) J,_i{n + N) 1 

J,{n + N + 1) q-^J,_i{n + N + 2) 



J^(n + 2A^ + 1) ■■■ J^(n + 3A^-1) q-^^ J^_i{n + 3N) 

Mn + 3) ■■■ J,{n + N + 1) q-^J,^i{n + N + 2) 

[-If ■: ... ; ; 

J^{n + 2N + l) ... J^(n + 3N-1) q-^^ J^^i{n + 3N) 
:-l)V' |3,,4,,...,iV + l„iV + 2;_i| , 



(46) 



which indicates that an additional factor (— l)^g ^ is needed. Taking this into account, 
we obtain 



rr+iH r^An + 3) - + 1) r^^in + 2) 

= -(l-g)«+iH^r'(^ + 3), 



(47) 



which is essentially the same as eq.(^). Similarly, eq.(^) is derived from the same identity 
as eq.(^) except that the column iV^_i is replaced by N^j^^. 
C. Equation (^) 



Finally we derive eq. (|2^) . We consider the identity of determinant, 

iu+i \ % 3,+i . . . 2N~~- i I 2iV + 01 











2N -\ 



v+l 



2A^+ 1 



(48) 

where ky, k^ and (pi are defined in eqs.(|3^, ( pTf ) and (|iT|), respectively. Applying the 
Laplace expansion to the right-hand side of eq.(H3), we have 



0= |2'„l+i,...,2iV^ W,2iV + l,+i| X |i,+i,l+i,...,2Aril,+i,0i| 
+ |i,+i,2;,3,+i,...,2iV^ Wl X |3,+i,...,2iV^ Wi,2iV + W,0i| 
- 3,+i, ■ . . , 2N^ 2iV + X |2'„ 3,+i, ■ ■ . , 2iV^ W, 0i| . (49) 



Using eqs.(P3), (BB) and (BPf), we get 



rj^^i(n + 2) rj^+^(n + 1) - q^'^^'^ {I + (1 - g)2g2"}^]^+iH ^^^(^ + 3) 
= (1 - g)g"+'^+^-^^(g'^ + q'^rNWin + 1) r;^(n + 2) . 



(50) 



Combining eqs.(pOD and (|20|), we obtain eq. (p^ . Thus we have proved that the r function 
(|I^) really satisfies the bilinear difference equations (PUD-(P^). 
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3.2 Derivation of dPni from Bilinear Difference Equations 

We now show that eq.([T8|) actually follows from the bilinear difference equations (|20D-( 
We first note that eq.(|l^) is factorized as 



"""^^^"-^ = {w^ + r-^n K - g-^-^^^) ~ ■ ^^^^ 

Let us rewrite the factors in eq.(0). Equation (|1^) gives 
which is rewritten by using eq. (^0[) as 
By the definition, it is obvious that 



Using eq.(E3D, we get 



From eqs.(|53| ) - (|55|) , we obtain 



-(1 - g)2g2n+2.+107V|^ _^ _ g)2g2(n-l)| 



/ A 



(53) 



(55) 



^ /^+i(n + l)r;^(ri + 2)r-+\(n - l)rj;(n) 

The factor q^^Wn — (1 — g)2g2n+i/+7Af jg rg-^yj^itten as follows, 
q'^^wn - (1 - g)^g^"+''+^^ 

^ r]^^,(n)rr(n + l) 
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where we have used eg .(1531) in the first hne and eg. (pi]) in the second hne, respectively. 



n+u+,N 'rN\\{n)TUn + 1) + (1 - q)q''-^''-^r'^+An)r''^^\n + 1) 



-U+l 

'N ' 



n 



Finally the factor q^'^Wn + (1 — q^q^"' '^+^^ ^ is rewritten as 

= (1 - q)q 
= (1 - 

X {(g-''+^-V;^+\(n - l)r;^(n + 2) - (1 - - l)rl 

+ (1 - {r'^,^,{n - lV^+\n + 2) + (1 - q)q--'r'^,X\{n - l)r^^{n + 2))} 

= (1 - {1 + (1 - ^f^''-''f^^^^^ , (58) 

where we have used eg. (^3]) in the first line and egs.(p2D and (p3D in the second line, 
respectively. From egs.(p7D and (|58|), we get 



- - (1 - (q'^'w^ + (1 - qf q'^^^''' 

= _(1 _ ^)2^2n+2.+10^|^ _^ _ g)2^2(n-l)| 

^ rj^+\(n + l)rMrNX\{n - l)r;^(n + 2) 

which is the same as eg. (^61) . Thus we have proved that the dPni eg.(pTD follows from the 
bilinear difference eguations (EH)- pi). 



4 Q'-Difference Analogue of Pni 

We have discussed the solution of dPni, and shown that it is expressed by the discrete 
Bessel function. It should be noted, however, that eg. (pT]) is essentially the same as the 
g-Bessel eguation, 

Mq^x) - {<f + q^")Ju{qx) + {1 + (1 - g) V} J,(x) = . (60) 



Actually, eg.(|TT]) yields eg.(^) by the replacement g" a;. This fact implies that the 
dPiii is a g-discrete system rather than a discrete system in essence. More precisely, let us 
consider the following eguation instead of the dPni (pTSD, 



w[qx)w[q x) = -—^ ^ , , . (61) 
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The solutions of eq. (|6l|) , with the parameters being the same as eq.(|T9|), are given by 



w{x) 



- Q 



(62) 



and the r function r^(x) is expressed as 



Ju{qx) 



(63) 



Here Jv{x) is the g-Bessel function of degree v. 

We next mention briefly the Lax pair for eq.(0). The auxiliary linear system for dPni 
is given by[|,|12[ 



<^n{qh) = Ln{h)^n{h) , 

<l>„+i(/i) = M„(/i)<l>„(/i) , 



(64) 
(65) 



whose compatibility condition is written as 



(66) 



If we choose suitable 4x4 matrices as L and M|T^, eq.(^) yields dPni (P^. The linear 
system for eq. (pT]) is written as g-difference system both in x and h, 



$(a;; qh) = L{x; h)^{x; h) , 
$(gx; h) = M{x; h)^{x; h), 



(67) 
(68) 



whose compatibility condition yields 



M(x; qh)L{x] h) = L{qx; h)M{x; h) . 



(69) 



Let the matrices L and M be 



L{x; h) 





w{x) 



\q'^w{q ^x) ' x'^ I h'^q^w{q ^x) 



h X 
c 



c 



wix) 



+ K 



C 



x^ 




wyx ) 
^2w{q~^x) 



X 

1^ 





A qwiq^^x) 

X 

w(x) , 



X 



A 



, (70) 
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M{x; h) 







/ {kx^ - q^\)w{x) 2 + c n \ 

7 — ^ 2 — X — 7 — r T U U 

q Aw{x) + cx q Aw{x) + cx 

10 

Hx'^ - X w{x) + /ix^ , (71) 
w{x) + A w{x) + A 

[ ^ yl 

respectively, where c, k, A and /i are arbitrary constants. Then we have from eg. (|69|) , 

w(qx)w(q~'x) = ^ (g'^^(^) + ^^')(^(^) + /^^') .72) 
"^^^"^^"^^^ q^{Kw{x) + c){w{x) + X) ' ^ 



which is essentially the same as eg. (]6T|) . 

Finally, we comment on the limit to Pni- We put w = {l—q)xu and take the limit q ^ 1. 
Then eg.(^) with the parameters (|19D reduces to the Pm eg.(0) with the parameters (|). 
Moreover, the r function ([63| ) reduces to that of Pni (0) up to a trivial factor and the 
dependent variable transformation ( |62D yields eg.(|D. 

In summary, 

1. Eguation (|6T|) reduces to the Pm eguation in the limit q ^ 1. 

2. Eguation (|61|) has particular solutions expressed by the determinant whose entries 
are given by the g-Bessel functions. 

3. The auxiliary linear system and the Lax eguation for eg.(pTD are written as q- 
difference eguations. 

These facts imply that eg.(|6lD is one candidate for the g-difference analogue of Pm. 

5 Concluding Remarks 

In this paper, we have investigated the solutions of special function type for the dPm 
eguation, and shown that they are expressed in terms of the Casorati determinants of the 
discrete analogue of the Bessel functions. The solutions are reduced to the Wronskian 
solutions of the Pm in the continuous limit. Hence we conclude that the dPm is a good 
discretization of the Pm. 

The discrete Bessel functions are transformed to the g-Bessel functions by a simple 
independent variable transformation. Based on this observation, we have presented a 
g-difference analogue of the Pm, whose particular solutions are expressed by the Wronski- 
type determinants of the g-Bessel functions. Moreover, we have mentioned that the Lax 
eguation and its auxiliary linear problem are written as g-difference eguations. 
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We expected that the r function of the dPni satisfies the bihnear form of the q- 



cyhndrical Toda molecule equation which we proposed in ref. since that of the Pni 
satisfies the cylindrical Toda molecule equation]]!^. However, we found it is not the case 
because of the asymmetry of structure of the r function. Actually the r function of the 
dPiii satisfies a slightly modified version of the g-cylindrical Toda molecule equation which 
also reduces to the cylindrical Toda molecule equation in the limit g — 1. Similar asym- 



metric structure of the r function has also been obtained for the solutions of the dPn[|T6 
We do not know yet whether such asymmetry is essential or not for the discrete Painleve 
equations. 

It is an interesting problem to study whether the other discrete Painleve equations 
proposed in refs. |^,[7|,p|, p!0| , pT| , p!2| admit similar determinant solutions or not. The existence 
of rational solutions together with their determinant structures is also an interesting subject 
to be investigated. 

Acknowledgment: We are grateful to Prof. Grammaticos and Prof. Ramani for fruitful 
discussions and useful comments. 
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